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We discuss a platform for the synthetic realization of key physical properties of helical Tomonaga 
Luttinger liquids (HTLLs) with ultracold fermionic atoms in one-dimensional optical lattices. The 
HTLL is a strongly correlated metallic state where spin polarization and propagation direction of 
the itinerant particles are locked to each other. We propose an unconventional one-dimensional 
Fermi-Hubbard model which, at quarter filling, resembles the HTLL in the long wavelength limit, 
as we demonstrate with a combination of analytical (bosonization) and numerical (density matrix 
renormalization group) methods. An experimentally feasible scheme is provided for the realization of 
this model with ultracold fermionic atoms in optical lattices. Finally, we discuss how the robustness 
of the HTLL against back-scattering and imperfections, well known from its realization at the 
edge of two-dimensional topological insulators, is reflected in the synthetic one-dimensional scenario 
proposed here. 


I. INTRODUCTION 

The rich interplay between orbital degrees of freedom, 
spin, and many-body correlations gives rise to fascinat¬ 
ing phenomena in quantum physics. New possibilities for 
their realization and observation are provided by the flex¬ 
ibility and control of quantum systems based on ultracold 
atoms in optical lattices m- Along these lines, syn¬ 
thetic magnetic fields HU or even effects of non-Abelian 
gauge helds such as spin orbit coupling (SOC) [7] have 
been observed in systems consisting of neutral atoms i- 
HI]. As we show below, this recent experimental progress 
even makes available the natural ingredients for the real¬ 
ization of quantum many-body systems which, while be¬ 
ing inspired by intriguing concepts from condensed mat¬ 
ter physics such as helical Tomonaga Luttinger liquids 
(HTLLs) [m [T3] , are not known to have a direct coun¬ 
terpart in real materials. 

In the context of strongly correlated one-dimensional 
(ID) Fermi gases, a new physical twist has been provided 
by the discovery of the HTLL mis]. In contrast to 
a conventional Tomonaga Luttinger liquid [TiHTT] . the 
HTLL is characterized by a peculiar locking of spin and 
direction of motion of the itinerant particles, namely that 
fermions with opposite spin move in opposite direction. 
The HTLL has been originally predicted as a metallic 
edge theory of two-dimensional topological insulators 
[IHHSI] exhibiting the quantum spin Hall (QSH) effect. 
Due to their robust spin sensitive transport properties, 
HTLLs are promising candidates for numerous spintron- 
ics applications. In our present work, we propose and 
study an exotic ID Fermi-Hubbard model that shows 
crucial aspects of HTLL physics at long distances, and is 
amenable to a systematic study and comparison to QSH 
edge states. The microscopic realization of this model in 
ID optical lattices draws intuition from engineered gauge 
fields HU- Reminiscent of SOC, our model contains a 



FIG. 1. a) Schematic of the two-component model Hamil¬ 
tonian § for b = an = 1, with hopping t between sites of 
the same component, and operations associated with a 

change in component and a hopping from left (right) to right 
(left), respectively, b) Schematic of model projected onto the 
lower band, resulting in an effective single component model. 
Grey ovals denote operators of the effective nearest neigh¬ 
bor hopping model in the lower band (c.f. Eq. @)- Physical 
lattice sites are denoted by dashed lines, c) Band structure of 
the Hamiltonian <§ for t = 6 = Qij = 1 and chemical poten¬ 
tial p at quarter filling. The coloring of the plots visualizes 
the momentum-dependent spin polarization of the fermions, 
where red denotes spin up and blue spin down. 


strongly spin-dependent hopping which, however, breaks 
time reversal symmetry, thus going conceptually beyond 
natural SOC. This unconventional feature is crucial for 
the basic physics characterizing the HTLL, namely that 
particles at opposite Fermi points have exactly opposite 
spin. 

The systems we investigate may be considered syn¬ 
thetic in a double sense. First, they are inspired by edge 
states of 2D QSH systems but are purely ID systems 
that do not rely on the presence of an insulating 2D 
bulk. Second, they are realized in synthetic material 
systems based on optical lattices in contrast to their 
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semiconductor based counterparts. The motivation of 
our study is hence also twofold. Regarding the first 
point, it is well known that a system which has no other 
degrees of freedom than a single helical channel cannot 
exist in a ID lattice system [^. We would therefore 
like to address both qualitatively and quantitatively the 
question as to which aspects of HTLL physics can be 
seen in a purely ID lattice system, and how much of 
the robustness against imperfections known from QSH 
systems survives in this synthetic scenario. Second, 
the experimental accessibility of subtle spin-dependent 
phenomena in HTLLs is limited in condensed matter 
systems by the presence of unavoidable imperfections, 
partly related to the surrounding bulk degrees of free¬ 
dom. By contrast, in synthetic systems based on cold 
atoms in optical lattices, the tunability of individual 
parameters and accessibility of observables with single 
site resolution may open up possibilities to measure, e.g., 
the spin correlation functions of the single channel HTLL. 


Below, we present an experimentally feasible scheme 
for the realization of our exotic ID Fermi Hubbard model 
with ^^^Yb atoms. The spin 1/2 degree of freedom of 
the particles in our model is encoded in two hyperfine 
levels with different magnetic quantum numbers mj. 
To engineer the hopping of the atoms, we use so called 
Raman assisted tunnelling techniques [IHS]. This allows 
us to imprint laser phases, and by means of dipole 
selection rules, also non-Abelian spin operations on the 
hopping processes. Beyond this experimental proposal, 
we extensively analyze our model theoretically. In the 
presence of on-site Hubbard interaction we assess, with 
a combination of analytical and numerical methods, in 
what sense its low-energy physics mimics a HTLL by 
computing long distance spin-spin correlation functions. 
Quite remarkably, the interacting model has a parameter 
line where these observable quantities can be exactly 
calculated analytically. Furthermore, we address the 
natural question to what extent the robustness against 
back-scattering - that can be understood as a topological 
protection in the natural realization of the HTLL as an 
edge state of a 2D topological insulator [T2J [TSl UHl 112] - 
carries over to our synthetic one-dimensional scenario. 


Outline - The remainder of this article is organized 
as follows. In Section |H| we discuss crucial features of 
the HTLL in general and introduce the microscopic lat¬ 
tice model that is at the heart of our present analysis. 
Thereafter, in Section [TTll we demonstrate how our model 
can be experimentally realized using state of the art tech¬ 
niques to control ultracold fermionic atoms in optical lat¬ 
tices. Section]^ is dedicated to the investigation of the 
low-energy physics of our model using a combination of 
analytical and numerical quantum many-body methods. 
Finally, in Section [V| we present some concluding re¬ 
marks and put our main findings into a broader context. 


II. MODEL BUILDING 

This section is concerned with the modelling of syn¬ 
thetic HTLLs. In Section [ll A[ we build a case for sim¬ 
ulating HTLL physics in ID lattice systems by summa¬ 
rizing the key observable features of the HTLL and com¬ 
paring the situation of its realization as an edge state of 
a 2D topological insulator to synthetic one-dimensional 
scenarios. Subsequently, in Section al B| we introduce the 
ID lattice model which forms the basis of our present 
study. 

A. Hallmarks of the HTLL 

The ideal HTLL has several intriguing properties 
which distinguish it decisively from both a spinless and a 
spinful single channel Tomonaga Luttinger liquid (TLL). 
Most prominently, there is no phase space to write down a 
spin-independent impurity term as discussed for the ordi¬ 
nary TLL in Refs. [231 mi- This is because of the absence 
of states with equal spin for particles moving in opposite 
directions. For edge states of the time reversal symme¬ 
try (TRS) preserving QSH state, this robustness can be 
pushed even further and precludes elastic single particle 
back-scattering due to any TRS preserving single particle 
term in the Hamiltonian I2niii2]- This makes the trans¬ 
port characteristics of the ideal HTLL quite unique and 
various subtle dissipative effects that can cause inelastic 
back-scattering even in the presence of TRS have been in¬ 
vestigated O [25H2^ . In a QSH system with a fixed spin 
quantization axis, the stability of the holographic HTLL 
at finite interaction strength has been investigated from 
first principles [50] . 

Another distinctive feature of the HTLL are its charac¬ 
teristic, very anisotropic spin-spin correlations. For spin 
operators S'-*- perpendicular to the spin quantization axis 
of the eigenstates at the Fermi momenta ±fcF, they read 
for large distances r as 

(S'^(r)S'^(O)) ~ cos(2fcFr)^, (1) 

i.e., they exhibit Friedel oscillations and follow a power 
law that depends on the correlation strength via the Lut¬ 
tinger parameter K. In contrast, for spin operators S" 
parallel to this quantization axis, they decay as 

(5"(r)5"(0))^l, (2) 

i.e., with an interaction independent power-law and 
without Friedel oscillations. Note that these correlation 
functions are a direct consequence of, and probe for, 
the helical nature of the Fermi surface: First, Eq. ([^ 
reflects that S" is a good quantum number around the 
individual Fermi points. This is because if S" were to 
mix left and right moving particles at the Fermi surface, 
these correlations would exhibit a non-universal power 
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law decay. Second, Eq. Q implies that spin-flips of 
S", as described by the S''*--operators, do couple the two 
Fermi-points, hence excluding that both Fermi points 
have the same polarization direction of S". 


Bloch Hamiltonian 

/lo(fc) = ^j. = 0,x,y,z 

d^{k) = (—tcos(A:),5cos(A:),ai^sin(/c),0). (4) 


In non-holographic, i.e., purely ID realizations of the 
HTLL (see, e.g.. Refs. piHHK] for various semiconductor 
based approximate realizations), there is a priori no 
topological protection against elastic single particle 
back-scattering by TRS and the extent to which the 
correlation functions concur with Eqs. has not 

been checked from first principles yet. In this work, we 
propose a minimal, experimentally feasible microscopic 
ID lattice model where the helical nature of the Fermi 
points is an exact feature at quarter filling. Regarding 
the characteristic spin-spin correlations, we are able to 
demonstrate the presence of a correlated HTLL state in 
the sense of Eqs. ([I][^ in a finite parameter range by a 
combined numerical (in the framework of density matrix 
renormalization group (DMRG)) [Ml I3Z] and analytical 
(combination of mapping to an exactly solvable model 
and bosonization) approach. Furthermore, turning to 
the mentioned robustness against back-scattering, we 
find that a spin-independent fluctuation in the lattice 
potential in our model gives r ise to a modified effective 
impurity term (see Section IV D| below for a more 


detailed discussion). This may be interpreted as a 
certain robustness of our synthetic HTLL which is hence 
also expected to have quite characteristic transport 
properties as compared to the conventional TLL the 
further exploration of which is an interesting future 
direction. We argue how the situation of an ideal HTLL 
can be more closely mimicked by slightly complicating 
the underlying band structure. Even though our model 
globally breaks TRS, there is an emergent TRS at or 
around the Fermi surface which gives rise to a certain 
robustness of transport properties. 


B. ID model for the synthetic HTLL 


The band structure and the Bloch functions explicitly 
read as 


E±{k) = ±\dl d = 


\u±{k)) = 


p±ikm 


(d^d^d^); 

= It) 


P±{k) = ^ (l±d(A:) -ct) 



(5) 


Due to the structure of d, we immediately see that at 
zLkp = ±|, i.e. at the Fermi points at quarter filling, 
the Bloch states of the lower band are characterized 
by (Ty| m_( i/ci?)) = ±|u_(±fcj’)), i.e., they are exact 
CTy eigenstates with opposite eigenvalue at the opposite 
Fermi points. This reflects the helical nature of the 
Fermi surface for a half filled lower band (quarter Filing 
of the lattice). In Fig. [^), we show the band structure 
of the free model for t — b — = 1 visualizing the 

lattice-momentum-dependent spin polarization of the 
Bloch states and indicating the helical Fermi points. 


We emphasize the role of the term bcos{k)ax appear¬ 
ing in the Bloch Hamiltonian Q. It differs crucially 
from an ordinary Zeeman term as induced by a magnetic 
field in a nano-wire (see, e.g.. Ref. [32]) because it is 
momentum-independent. A Zeeman term Bax opens a 
gap between the two spin species, even at /c = 0 where 
the Rashba term siii{k)ay vanishes. However, it also 
gives an equal ax spin polarization to the Bloch states at 
opposite Fermi points thus competing with the desired 
helical nature of the Fermi surface. In contrast, the term 
bcos{k)ax appearing in Eq. Q is maximal at fc = 0 but 
vanishes exactly at k = ±|, i.e., at the Fermi points at 
quarter Filing. The Rashba term being an odd function 
of momentum then gives rise to exactly helical Fermi 
points. 


We consider a lattice model with a single spin- 
i fermionic degree of freedom per site and choose the 
lattice constant as our unit of length. The free tight- 
binding model of the underlying band structure is given 
by 

Po = [b<Jx + ioiRCTy - ta^j] '0J+1 -I- h.c., (3) 

3 

where tpj = (V'i.T’spinors of fermionic Feld 
operators, t is the ordinary spin-independent hopping 
strength, an denotes the Rashba velocity, and b tunes 
a TRS breaking hopping term which may be seen as an 
exotic Zeeman term. On Fourier transform we obtain the 


We note that bcos{k)ax could be replaced by 
6 cos"(fc)tTa; in Eq. Q at the expense of introducing 
n-nearest neighbor hopping in Eq. (|^. That way, 
not only the Fermi points but also an expansion to 
order around them would be exactly helical. For 

n = 2, for example, the full linearized theory around the 
Fermi surface would give an exact HTLL up to band 
curvature terms. In the following, we focus on the case 
n = 1, where the model is to some extent amenable to 
analytical study even in the presence of interaction. 

Due to the absence of a stable Fermi liquid in ID, con¬ 
sidering interaction effects is of key importance to make 
experimentally relevant predictions on helical liquids. 
Here, we model the interaction between the fermionic 
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atoms with an ordinary on site Hubbard interaction term, 

Hi = U^ '^3,'y = cr =t, i (6) 

j 

with interaction strength U. In Section [TV] we show that 
the model described by 


H = Ho + Hi (7) 

indeed exhibits HTLL physics in a finite parameter range. 


III. IMPLEMENTATION WITH COLD ATOMS 
IN OPTICAL LATTICES 


Intuition for the practical realization of our model ^ 
is drawn from both theoretical proposals 00 and recent 
experiments on synthetic classical gauge fields 

in optical lattices, where Raman assisted tunnelling 
techniques are used to engineer the phase of the hopping 
amplitude of the atoms. The experimental scheme 
proposed here is related to the experimental setup in 
Ref. m, where a magnetic flux in a so called synthetic 
dimension formed by internal states of ^^^Yb atoms 
has been realized. However, as we detail below, our 
present proposal goes crucially beyond this scheme since 
it contains a spin-flip hopping process, which would 
correspond to a combined hopping in the synthetic 
dimension (i.e. the spin flip), and hopping in the real 
dimension (see diagonal hopping terms in Fig. [^a). 


We now detail how the tight-binding model (§ can 
be realized using Raman assisted tunnelling techniques 
in a system of ultracold fermionic atoms in a ID optical 
lattice. Motivated by recent experimental progress m, 
we focus on the fermionic alkaline earth atoms ^^^Yb. 
The ground state in this scenario is a state, with 
F" = / = I for ^^^Yb, and we choose to form the 
spin a in Eq. (|^ with mj = —5/2 and m'j = —1/2. 
The occupation of the remaining states in the manifold 
rrif = —5/2,..., -1-5/2 is inhibited by a combination of 
dipole selection rules and energy conservation arguments, 
as we will explain in more detail below. This effective 
projection onto the subspace forming a has already been 
experimentally achieved in Ref. m- We note that the 
most relevant parameter regime for our proposal is de¬ 
fined by 6 = an in Eq. <§ , where our model is also 
amenable to analytical study (see Section IV). However, 
our implementation will also allow for the case b ^ an- 
We start by lifting the degeneracy of the six TO/-states 
by a Zeeman splitting Aj, induced via a magnetic field 
in the z-direction. Then, to implement the Hamiltonian 
Eq. (§ for h = an^ we need to engineer the following two 
nearest neighbor hopping terms (see Fig. 


I. A spin-dependent hopping of strength b = an, 
where hopping from left to right is accompanied 
by the spin operation S+ = ^{<Tx + icTy) while the 


E 



- 1 / 2 ) 

-5/2) 


FIG. 2. Proposed setup for the implementation of the spin- 
dependent hopping in Eq. ([^ at the point an = h using Yb. 
The two spin states are given by the nuclear zeeman states 
of the ^So ground state, |m/ = — |),|m/ = — |), subject to 
a Zeeman splitting of 2Az. Raman assisted tunneling with 
lasers with Rabi frequencies f2i/2, frequencies uJ- 1/2 and polar¬ 
izations respectively, couple nearest neighbor sites via 

a I^Pi) virtual excited state, detuned by an amount 5. The 
lattice is tilted with equal tilting A of the two spin states. 


hermitian conjugate S- acts when hopping from 
right to left. 

2. A spin-independent nearest neighbor hopping of 
strength t. 

We start by addressing the first term. This term is 
related to what has been demonstrated in m, however 
again we emphasize the difference here is that the spin 
flip term (hopping in the synthetic dimension) is associ¬ 
ated with nearest neighbour hopping (in the spatial di¬ 
mension). This additional ingredient is physically crucial 
here to achieve the exact helical polarization of the Fermi 
points in our model and is hence at the heart of our cur¬ 
rent proposal. To implement this specific nearest neigh¬ 
bor hopping term, the natural spin-independent hopping 
J stemming from the overlap of the Wannier functions 
is prohibited by tilting the optical lattice with a slope 
A ^ J. Several ways to achieve this have been reported 
0. Here, we need to tilt the two spin species with an 
equal slope, which speaks against using a magnetic field 
gradient. Instead, ways of achieving this are to employ a 
spatial gradient in the laser intensity that may be real¬ 
ized by a laser focus, to ’’shake” the lattice by means of a 
time-dependent frequency difference of the laser beams, 
or to simply tilt the lattice by means of the gravitational 
potential. This can generate a slope in the optical lattice 
potential which is equal for our two spin species. 

The hopping is then restored in the tilted lattice with 
the help of two far detuned Raman lasers [5]. The locking 
between the spatial hopping direction and the spin flip 
operators S±, respectively, is achieved by dipole selec¬ 
tion rules. More concretely, due to nrifi — ruf = 2h, two 
units of angular momentum along the z-direction need 
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to be transferred in every hopping process. To this end, 
we use a pair of Raman lasers with frequencies uji,uj 2 , 
polarisations and wave vectors CTi = tT_)_, ki = and 
CT 2 = CT-, k 2 = respectively. The Raman lasers 

are detuned by 6 from a manifold of excited states, here 
given by the state of ^’^^Yb. For the hyper- 

fine splitting between the excited states with different 
F' = ( 5 , 1 ) is about 4.7 GHz [35], so that the detun¬ 
ing S of the Raman process may also be on the order of 
1 GHz [ID, i.e., far detuned compared to both the en¬ 
ergy scale of the optical lattice and the line-width of the 
corresponding optical transition. If the frequency differ¬ 
ence concurs with the tilting and the Zeeman splitting as 
UJ 1 —UJ 2 = A-I- 2 A 2 , spin-flip nearest neighbor hopping as¬ 
sisted by a two photon Raman process conserves energy. 
As well, hopping processes corresponding to a climb (de¬ 
scend) in the tilted lattice are associated with distinct 
photon processes, namely oji is absorbed (emitted) and 
UJ 2 is emitted (absorbed), respectively. Due to the dis¬ 
tinct polarizations cti = a+,a 2 = ct- of the photons, 
a climb (descend) in the lattice corresponds to a trans¬ 
fer in the angular momentum along the z-direction of 2h 
(—2h), precisely matching the difference in the magnetic 
quantum numbers between our two spin species. This 
establishes a one-to-one correspondence between the di¬ 
rection of the Raman assisted hopping and the applica¬ 
tion of the spin flip operators S± (see Fig. and pre¬ 
vents the dynamical occupation of three of the other four 
m/-states m/ = —3/2,1/2, 5/2 that are not part of our 
spin (T. More concretely, rrif = —3/2,1/2,5/2 are en¬ 
ergetically off-resonant by at least A^ and, in addition, 
cannot be reached from our spin states by transferring 
an even number of units of angular momentum along the 
z-direction. In order to prevent occupation of the fourth 
state, m/ = 3/2, its energy is shifted away from reso¬ 
nance by a light shift El. 

In order to not occupy higher Bloch bands of the op¬ 
tical lattice, the hopping strength aji resulting from our 
Raman assisted spin-flip tunnelling scheme should be 
smaller than band gap Af, between the lowest and the 
second lowest Bloch band. Lastly, we note that unwanted 
onsite spin-flips, which would correspond to an ordinary, 
i.e., momentum-independent Zeeman term in (§ are off- 
resonant by the lattice tilting A and are hence strongly 
suppressed if aji < A. 

We point out that the model engineered so far, i.e., 
Eq. ([^ with f = 0 at half filling realizes a flat band ID 
topological insulator, similar to the model introduced by 
Su, Schrieffer, and Heeger [331. However, in order to ob¬ 
tain a finite Fermi velocity at the helical Fermi points at 
quarter filling, a finite t is necessary to make the flat 
bands dispersive. This brings us to the second term 
above, the spin-independent hopping t. This can be re¬ 
alized by another two-photon Raman assisted tunnnel- 
ing process which is not spin-selective. However, for this 
process, a spin-dependent phase may be implemented by 
varying the direction of the Raman lasers in order to give 
their wave vector a component parallel to the ID optical 


lattice (^-direction). This allows us to tune away from 
an = b, thus enabling the implementation of Eq. ([^ in 
all generality. 

Finally, we note that the helical nature of the Fermi 
points occurs at quarter filling of the lattice, i.e., at 
a particle density of one atom per two lattice sites on 
average. This commensurate filling could be achieved 
experimentally by, for example, temporarily switching 
on a superlattice with the double period and realizing a 
Mott phase with one particle per site in this super lattice. 


IV. ANALYTICAL AND NUMERICAL 
ANALYSIS 

Our Hamiltonian Q represents a ID Fermi-Hubbard 
model with a somewhat exotic underlying band struc¬ 
ture (at & = Ofl = 0 it would be the ordinary Hubbard 
model). The long wavelength physics of this model will 
be revealed with a combination of analytical and numer¬ 
ical methods in the following. 


A. Exact solution and lower band projection 


As a first approach, we present an elegant way to map 
the low-energy physics around the Fermi points of this 
model to that of an exactly solvable model. We con¬ 
sider the lattice at quarter filling and put b = aii,t = I. 
We project the model to its lower band which gives a 
good approximation of the low-energy physics around the 
Fermi surface, if the energy scales of interest are much 
smaller than the energetic separation 2& — 1 of the Fermi 
surface from the upper band (see Fig. [fright panel). 
The effective free Hamiltonian then reads as 

( 8 ) 

k 


where are the creation operators of the Bloch states 
in the lower band, i.e., Zj(|0) = |u_(fc)). Now, we also 
project the interaction term 77/ to the lowest band. To 
this end it is helpful to have a Wannier basis of localized 
states spanning the lower band. The /c-independence of 
\d\ at b = aji significantly helps here (see Eq. (§)). Fol¬ 
lowing Ref. [40] and considering that the ordinary hop¬ 
ping term t does not influence the Bloch states due to its 
spin-independence, we obtain 


= +I1.C. 

3 


(9) 


where we have shifted the Energy by b units to get 
rid of the constant term \d{k)\ = b in the spectrum 
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E-{k) = — cos(A:) — \d{k)\ of the lower band. The effec¬ 
tive momentum-dependent model for the lower band 
is visualized in Fig. [^). The projection Hj of the in¬ 
teraction Hamiltonian H: to the lower band assumes in 
terms of the localized Wannier operators 7 ^ the simple 
form 


U 


3 

Putting together Eq. (TpOl, we have 


( 10 ) 


H =V-HV- =Ho + Hi = 


XI [t'It'I+i + ^ fijhj+x, (11) 


where V- denotes the projection onto the lower band. 
Eq. (11) represents an effective spinless fermion model 


in terms of the operators 7 ^ with only nearest neighbor 
interaction which is commonly referred to as the t — V 
model. However, via the spin structure of the 7 ^ 
operators (see Eq. a spin dependence is encoded in 
this effective spinless model. Here, t = 1 and 17/4 plays 
the role of V. This model can be solved exactly using 
the Bethe ansatz and is, at half filling, known to describe 
a Tomonaga Luttinger liquid (TLL) for interactions 
U < A. At [/ > 4 the system develops a 2kp charge 
density wave order and a gap is opened. Note that the 
quarter filling of the full lattice we started with now 
corresponds to half filling of the effective one-band model. 


B. Bosonization and characteristic spin correlations 

Our goal is to verify that the long distance physics of 
our model indeed exhibits the characteristic behaviour of 
the HTLL. The most striking signature of the HTLL is 
the asymptotic decay of the spin-spin correlation func¬ 
tions which depends drastically on the spin direction. 
The spin-spin correlations perpendicular to the polar¬ 
ization axis of the eigenstates (here the j/-axis) decay 
with a non-universal power law that depends on the Lut¬ 
tinger parameter K and exhibit Friedel oscillations (see 
Eq. (§). In contrast, the spin-spin correlations parallel 
to this axis decay with the second power in distance, in¬ 
dependent of the interaction strength, and do not show 
Friedel oscillations (Eq. (§). From Eq. ([^ we know that 
the |u_(±fci?)) at Ekp = ±| describe perfectly helical 
modes indicating that the long wavelength physics of our 
model resembles the HTLL. However, in the presence of 
interaction, long range correlation functions of the pro¬ 
jected model in Eq. (10) are hard to access from its exact 


Bethe ansatz solution. Still, we know from the exact so¬ 
lution that the low-energy theory of the effective spinless 
model (11), is a spinless TLL the correlation functions 


of which are analytically computable. We hence project 
the microscopic lattice spin operators Sj, i = x,y,z at 


site j to the lower band and bosonize them in a long 
wavelength continuum model that is linearized around 
the Eermi energy. To this end, we first decompose the 
Eermi operators 7 j of the lower band into a right-moving 
and a left moving part 


Ij = iPl{.x) + iPr{x), x=j 


( 12 ) 


and treat x as continuous parameter in the subsequent 
analysis. Following the notation of Ref. [H], the 
bosonized form of the operators reads in the thermody¬ 
namic limit as 


_ ipkpx -i{p4>{x)-B(x)) 


^pix) = 

y Zira 


, p = R/L = ±, 
(13) 


where Up are the mutually anti-commuting Klein factors, 
a is a short distance cutoff, and </, 9 are the bosonic phase 
field and its dual, respectively. The main merit of the 
bosonized representation is that the interacting model 
(11), linearized around the Fermi energy, is quadratic in 


the bosonic fields which allows for the analytical calcu¬ 
lation of the long range correlation functions (see, e.g.. 
Ref. [41]). To employ this for the calculation of the 
desired spin-spin correlation functions, we project the 

lattice spin operators S'® = i = x,y,z to 

the lower band. Looking at Eq. ([^, it is clear that this 
amounts to the mapping 






(14) 


Along with Eqs. (12), (131, the spin operators can now in 
principle be brought into bosonized form. However, from 
Eq. (14), it is clear that the off-diagonal spin operators 
Sj, i = x,y, contain lowest band operators 7 j, 7 j_i at 
neighboring sites. Explicitly, 

Sj ^ -^(tX-i +7]-i7j), ^ - ih-i)- 

In the linearized continuum model, by virtue of 


Eq. (12), operators like 'ip^ix — l)ipp'{x), p,p' = L,R 
hence appear which are separated by the short distance 
of 1 < a. In order to evaluate the relevant correlation 
functions using the bosonic long wavelength theory such 
terms first need to be simplified to operators evaluated 
at a single position by performing an operator product 
expansion (OPE). Terms with p = p' above give lead¬ 
ing contributions proportional to or d^d. Pair cor¬ 
relations evaluated at positions x and y of such terms 
give rise to the universal power-law \x — ?/|“^ appear¬ 
ing in Eq. (|^. In contrast, the leading OPE of terms 
with p ^ p', e.g., ip^pfix — 1 )iPl{x) gives rise to operators 
that oscillate with 2kFX and that depend on the bosonic 
fields as _ Pair correlations of such operators give 

rise to the non-universal power-law \x — y\~‘^^ appearing 
in Eq. 0 and their oscillatory behavior causes the con¬ 
comitant Friedel oscillations. The explicit evaluation of 
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the long distance spin-spin correlators {S^{x)S^{y)) with 
~ y\ ^ 1; is now tedious but straight forward. For 
the parallel {S^{x)S'^(y)) correlations, relative signs in 
the matrix structure of the spin operators lead to a can¬ 
cellation of all oscillatory terms, hence resulting indeed in 
the universal power law and the absence of Friedel oscil¬ 
lations displayed in Eq. (H). For the perpendicular spin- 
spin correlations, in contrast, the leading correlations at 
repulsive interactions {K < 1) are produced by the os¬ 
cillatory terms which cause the 2k f Friedel oscillations 
and give rise to the non-universal power-law appearing in 
Eq. 0. The bosonization analysis hence confirms that 
the spin-spin correlation functions of our model at quar¬ 
ter filling and b = a ^ t decay with the characteristic 
behavior of the HTLL. At iF = 1/2 two particle umklapp 
scattering terms become relevant and the system enters 
a charge density wave phase as predicted by the Bethe 
ansatz solution. 


C. Numerical analysis 


Complementary to the previous analytical approach, 
we treat the full microscopic lattice model Q in the 
framework of finite system DMRG m- This allows 
for a direct measurement of the correlation functions 
along the directions parallel and perpendicular to the 
polarization at the Fermi surface. More specifically, 
we measure various spin-spin correlation functions at 
positions ri and r 2 and fit their decay of to the power-law 
with r = |ri — r 2 |. To minimize finite size effects 
we keep ri and r 2 at fixed fractions of the system length 
L, namely ri = 3L/4, r 2 = Lj'^ and vary the distance 
r = L/2 by variation of the system size L [30]. Besides 
confirming our analytical analysis from first principles, 
our numerical study also allows us to access a broader 
parameter range where the model is not amenable to 
analytical treatment. 

We first concentrate again on the case b = an. For 
the parallel spin direction, we find that the correlations 
are basically independent of interaction, confirming the 
behavior given in Eq. 0. By contrast, the correlation 
functions in the perpendicular direction depend heavily 
on the interaction strength, in agreement with Eq. Q. 
The dependence of the Luttinger parameter K on the 
interaction strength U can be extracted from finite size 
scaling of the perpendicular {S^{ri)S^{r 2 )) correlations, 
where K = ^. Numerical errors from our DMRG 
calculations and finite size effects limit the precision of 
the fits for ^ to a precision of a few percent. The results 
are shown in Eig. for separation of bands increasing 
from b = an = t to b = an = 50t. When the band 
separation is much larger than the interaction strength, 
the projection to the lower band is a good approxima¬ 
tion, and the Luttinger parameter follows the behavior 
of the Bethe Ansatz [42] within the numerical errors of 
our finite size scaling. In this limit, the system crosses 
a phase transition into a charge-density wave phase at 



FIG. 3. The decay exponent ^ extracted from a fit to 
|ri — r 2 |~^^ as a function of the Hubbard interaction U, as ex¬ 
tracted from DMRG measurements of the correlation function 
{S^(ri)S^(r 2 )) (solid) and {S^(ri)5'^(r2)) (dashed), respec¬ 
tively. Here, ri = 3Z//4 and r 2 = L/A, where the system size 
L is varied between 40 and 140 sites. Different colors corre¬ 
spond to separation of the bands increasing from b = an = t 
to b = an = 50t. Additionally, the solution of the Bethe 
Ansatz is shown for comparison. The value of ^=0.5 indi¬ 
cating the phase transition between the Luttinger liquid and 
charge-density wave phase is shown with a dashed black line 
as a guide for the eye. The spread around ^ = 1 at (7 = 0 
gives an indication of the magnitude of the numerical error. 
Inset: Gorrelation function (S'^(ri)S'^(r 2 )) as a function of 
\ri — r 2 \, and the best fit from which the exponent ^ can be 
extracted. Here all data is at an — b = lOt, and the colors 
indicating U — 0 increasing to 17 = 4t. 


17 = 4; the critical point is at a Luttinger parameter of 
K = 0.5 |33|. However, when the interaction becomes 
comparable to the band separation, projecting out the 
upper band is no longer a good approximation as the 
interaction can readily mix population in both bands, 
and the behavior of the system deviates systematically 
from the prediction of the Bethe Ansatz. 

We now consider the case b ^ an, where the disper¬ 
sion of the lower band is not just a simple —fcos(fc) but 
is influenced by the spin-dependent hopping terms (see 
Eq. if)- As a consequence, the lowest band projection 
is no longer given by the exactly solvable model ( [TT] ). 
However, as long as \b — an] <t, there are still only two 
exactly helical Eermi points at quarter filling. In the up¬ 
per panel of Eig. [^ we show the [/-dependence in the 
decay exponent ^ of the parallel and perpendicular spin- 
spin correlation functions |ri — r 2 |~^^ for fixed an = 5t 
and various values of b. The concomitant dispersion of 
the lower band for each value of b is shown in the lower 
panel of Fig. [^ The qualitative structure of the correla¬ 
tion functions is still similar to the case b ^ an, in par¬ 
ticular the striking anisotropy between the parallel and 
the perpendicular correlations. Discrepant points in the 
exponent of the perpendicular correlation function ap¬ 
pear when the difference an — b becomes relevant when 
compared to the interaction strength. However, when 
the interaction is strong enough with respect to an — b, 
the behavior is similar to that of the case an = b. We 
hence conclude that the HTLL long distance behavior of 












1.25 


be microscopically modeled by the operator 



FIG. 4. Upper Panel: Decay exponent ^ from fit to |ri —r 2 |“^^ 
as a function of the Hubbard interaction U, as extracted 
from DMRG measurements of the correlation function paral¬ 
lel {5'"(ri)S"'(r2)) (dashed) and perpendicular (5^(ri)S'^(r2)) 
(solid) to the spin quantization axis. Here, ri = 3I//4 and 
r 2 = D/4, where the system size L is varied between 40 and 
140 sites. Different colors correspond to 6 = 4.2t to 6 = 5t, 
with Qij = 5t fixed for every value of b. Lower Panel: Band 
structure of the lower band for the parameters of the upper 
panel. 


our model persists also in an extended parameter regime 
where it is not exactly solvable. 


D. Robustness against back-scattering 

A hallmark of the HTLL is its robustness against 
single particle back-scattering. For the ideal HTLL con¬ 
sisting of nothing but a right moving branch with spin 
up and a left moving branch with spin down conjugated 
by TRS, elastic single particle back-scattering cannot 
be induced by any TRS preserving perturbation. In 
particular, a generic spin-independent impurity term 
causing back-scattering in an ordinary TLL cannot exist 
since the spin needs to be flipped in order to couple 
the opposite branches. In contrast, our present model 
([^ is a microscopic lattice model with a spin degree of 
freedom for which any local operator, irrespective of its 
spin dependence can be written down. The question that 
we address in this Section is to what extent phenomena 
reminiscent of this protection against back-scattering 
occur in our synthetic model, which, away from the 
Fermi surface, deviates from the ideal HTLL. 

The simplest conceivable imperfection, a spin- 
independent impurity with strength Vq at site j = 0 can 


Oimp(j = 0) = VbV’JV'o = Vb(V'(5,^V'o,t + (15) 

The Fermi points at zLkp = ±f have exactly opposite 
spins and thus the matrix element for scattering from kp 
to —kp vanishes. The admixture of the opposite spin 
polarization away from the Fermi surface scales linear 
with the deviation q = k — kp. To assess the effect of 
Oimp on the low-energy theory around the Fermi surface, 
we first perform a projection to the lowest band. With 
= Ijlj and ^ (see Eq. 0), we 

obtain 


Ou 




(16) 


Using Eq. (12), we translate Oimp into the long wave¬ 
length continuum theory, which gives 


u- 


Dlimp — Oimp(j) — ^ + hj_l) 


(17) 


dxV(x) p(x) + (e -I-h.c.^ 


where we defined p{x) = 'tjj\^{x)'tjjp{x) + tjj]^{x)ipfi{x) and 
kp = ^ has been used. In bosonized form, the potential 
term Himp reads as 



xV{x) 



1 

27ra 



i(2kpx-2(l>) ^ ^ 


where the first term describes forward-scattering and the 
second term back-scattering. Compared to the conven¬ 
tional spinless TLL, where the back-scattering of a po¬ 
tential is described by a term ~ V(x) cos(20(a::)) P51IM] . 
there is an additional pre-factor of dx4>{x) in our model. 
Simple power counting would hence indicate that back- 
scattering is less relevant in our model. However, to 
higher order in perturbation theory, the combination of 
forward-scattering and back-scattering may also generate 
the conventional term ~ cos{2(j){x)) in our model, but 
with leading order in the potential strength. Thus, in 
the perturbative regime, where temperature and/or bias 
voltage are larger than the potential strength V, back- 
scattering is expected to be suppressed in our model com¬ 
pared to a conventional spinless TLL. 

We note that if we modify our model (§ by sub stitut - 
ing bcos^{k) for bcos{k) (see discussion in Section HB|, 
the impurity operator given in Eq. ( [T^ becomes less 
relevant in the effective many-body theory and enters 
only via band curvature terms which go beyond the 
linearized TLL model. 


We stress that the above discussion assumes that the 
Eermi wave vector is tuned exactly to = ^. In prac¬ 
tice, the Fermi energy may slightly deviate from this 
value, say kp—^ = 6. This will give rise to a conventional 
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impurity term ~ 5V{x) cos(2(^(a;)) on the order of this de¬ 
viation. However, if the bias voltage Ub ot the tempera¬ 
ture T are larger than 5V, this term will not pinch off the 
conductance. In conclusion, for < Ub,T <V, a 

scattering potential of strength V would pinch off a con¬ 
ventional TLL while the synthetic HTLL is only pertur- 
batively affected. In this sense, the experimental verifi¬ 
cation of the reduced back-scattering discussed here does 
not require fine tuning of the Fermi momentum. 

V. CONCLUDING REMARKS 

We have demonstrated how HTLL physics naturally 
emerges in a synthetic framework based on ultracold 
fermionic atoms in optical lattices. In this context, static 
spin-spin correlation functions have been calculated as a 
hallmark of HTLL physics. These observable quantities 
are experimentally accessible in state of the art exper¬ 
iments on optical lattices and provide a probe of the 
helical nature of the Fermi surface. Furthermore, the 
robustness of the ideal HTLL against back-scattering 
has been shown to have an interesting counterpart in our 
synthetic realization. In a broader context, the study 
of mesoscopic transport properties of neutral atoms in 
optical lattices is an emerging field of research (see, 
e.g., Ref. mi)- Along these lines, our present work 
may open up a playground for the study of intriguing 


transport properties in the framework of unconventional 
quantum impurity problems. Once the single channel 
regime becomes accessible in transport experiments on 
ultracold atoms in optical lattices, HTLL physics in 
the proposed setup may also be dynamically probed by 
observing the correlation between the spin polarization 
and propagation direction of excitations. 
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been reported in a model with spin-independent hopping 
which maps to a Rashba spin orbit coupled wire. 
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